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We have developed a macroscopic description of coherent electro-magnetic radiation from air 
showers initiated by ultra-high energy cosmic rays due to the presence of the geo-magnetic field. 
This description offers a simple and direct insight in the relation between the properties of the air 
shower and the time-structure of the radio pulse. 
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I. INTRODUCTION 

In recent years the interest in the use of radio detection 
for cosmic ray air showers is increasing with the promis- 
ing results obtained from recent LOPES [l|, 0| and CO- 
DALEMA Q experiments. These experiments have in 
turn triggered plans to install an extensive array of radio 
detectors at the Pierre Auger Observatory [4|. There is 
thus a growing interest the link between the properties 
of the air shower and the time structure of the emitted 
pulse. Already in the earliest works on radio emission 
from air showers [!, 0, @, [U, the importance of coher- 
ent emission was stressed. Two mechanisms, Chcrenkov 
radiation and geo-magnetic radiation were proposed as 
possibilities. In more recent work 

[HGUm, the picture 
of coherent synchrotron radiation from secondary shower 
electrons and positrons gyrating in the Earth's magnetic 
field was proposed. Extensive results on geo-synchrotron 
emission, based on realistic Monte-Carlo simulations of 
the shower development, are given in (HI. 

The primary motivation of this work is to improve on 
the understanding of the relation between the measured 
pulse shape using radio receivers and the properties of 
the air shower induced by a cosmic ray. Therefore, we 
performed macroscopic calculations which allow, under 
simplifying conditions, to obtain a simple analytic ex- 
pression for the pulse shape. This analytic expression 
shows a clear relation between the pulse shape and the 
shower profile. 

The picture we use is very similar to that used in 
Ref. Q which we refine by using a more realistic shower 
profile and where we calculate the time-dependence of the 
pulse. The magnetic field of the Earth induces, by pulling 
with the Lorentz force the electrons and positrons in op- 
posite directions, a net electric current in the electron- 
positron plasma. This plasma moves with almost the 
velocity of light towards the Earth at the front end of 
the cosmic-ray air shower. In our approach the collec- 
tive aspect is emphasized by treating this induced cur- 
rent as macroscopic. This differs from the approach of 
Refs. [TU, [l3[ where the motion of individual particles 
is stressed (microscopic approach). In both the macro- 
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scopic and the microscopic approach the emission of the 
electromagnetic pulse is caused by moving charges in 
the Earth's magnetic field. Therefore these two pictures 
should be regarded as presenting a complementary view 
of the same physical phenomenon. There are however 
differences in the predicted pulse shapes and we hope 
that by presenting this complementary picture the un- 
derstanding of radio emission from extensive air showers 
can be improved. 

In the introduction of Section |TT] the basic outline of 
our approach is presented and the various aspects are de- 
tailed in the different subsections. Starting from a very 
basic picture we present our results in Section Hill Subse- 
quently the effects on the pulse shape are investigated of 
finite lateral extend, finite pancake thickness, and a re- 
alistic energy distribution of the electrons and positrons 
in the air shower. 



II. THE FORMALISM 

When an UHE cosmic-ray particle enters the upper 
layers of the atmosphere, a cascade of high-energy par- 
ticles - called a cosmic-ray air shower - develops. Due 
to the high velocities, most of the particles are concen- 
trated in the relatively thin shower front, which, for obvi- 
ous reasons, is called the 'pancake'. The pancake, which 
for the present discussion is assumed to be charge neu- 
tral, contains extremely large numbers of electrons and 
positrons. Near the core of the shower this pancake has 
a typical thickness of a few meters and is moving to the 
surface of the Earth with (almost) the velocity of light 
through the magnetic field of the Earth. The Lorentz 
force on the charged particles induces an acceleration of 
the particles in the x direction, which is perpendicular to 
the magnetic field and the shower axis. However, due to 
the frequent collision with the air-molecules, where the 
relatively small transverse velocity is randomized, this ac- 
celeration, when averaged over all electrons, rather trans- 
lates into a drift velocity and thus an electric current in 
the x direction. This picture is similar to what happens 
to electrons in a copper wire. When a voltage is applied 
over the wire, the electrons undergo a constant accel- 
eration due to the electric force which is however com- 
pensated by collisions with the copper atoms, resulting 
in a constant drift velocity and thus a constant electric 
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current. At the surface of the Earth, electromagnetic ra- 
diation can be detected, which is due to this relatively 
constant electric current moving with high velocity to- 
wards the Earth. The shape of the electromagnetic pulse 
is principally determined by the (relatively slow) varia- 
tion in time of the magnitude of the current, combined 
with time retardation effects. 

To emphasize the basic principles, we confine ourselves 
to a rather simple geometry where the cosmic shower 
moves straight towards the Earth's surface (the — z direc- 
tion, see Fig. [T]) with velocity v s = —f3 s cz where (3 3 « 1. 
The position of the shower front above the Earth's sur- 
face is given by z = — (3 s ct, where the front of the shower 
reaches Earth at time t = 0. The Earth's magnetic field 
(with magnitude Be) is parallel to the surface (in the 
y direction), B = Be])- The strength of the induced 
electric current depends on the distance h from the front 
of the shower and on the time t in the shower develop- 
ment. The direction of the current is in the x direction. 
All quantities are measured in the rest system of the ob- 
server who is at rest at the surface of the Earth. 
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FIG. 1: [color online] The basic geometry as used in this paper. 
The pancake is shown as the shaded area, moving with velocity 
v B . The observer is at a distance d from the point of impact. 
The curvature of the pancake has not been included in the 
present calculations. 

In the calculation of the current density we will ini- 
tially assume a finite extent in the horizontal directions 
(x and y). However, soon we will integrate over these 
variables, knowing that the charged-particle density is 
strongly peaked near the center of the shower. To empha- 
size the importance of the distance h behind the shower 
front, we will write the electron/positron density as 



p e (x, y, z,t) = J p e (x, y, z, t, h)dh , (1) 

where we assume a simple factorized form, 
p e (x,y,z,t,h) = N e d(z + /3 s ct + h) 



x ft(i)pNK.a(x,y)p p (h) 



(2) 



The total number of charged particles at the time of max- 
imum shower development is denoted as N e . The veloc- 
ity of the shower front is given by v s — ~P s cz. The 



lateral distribution function is normalized according to 
/ Pnkg(x , y) dx dy = 1, the pancake distribution obeys 
a similar normalization, J p p (h) dh = 1, and the maxi- 
mum of the temporal (or longitudinal) distribution f t (t) 
is normalized to unity. A detailed discussion of the pa- 
rameterizations for these shower functions is given in the 
appendix. In Section IA 41 also the effects of an energy 
spread of the electrons and positrons are considered. 

To emphasize the collective aspects of the model the 
calculation of the drift velocity of the electrons and 
positrons is treated as a separate topic. The magnitude 
of the induced current is calculated as the number of 
electrons (and positrons) multiplied by an average drift 
velocity. In the following stage this is combined with the 
shower profile to calculate the emitted electromagnetic 
pulse. 



A. Magnitude of the Current 

For the present estimate it is assumed that there are 
equal numbers of positive and negative charges moving 
towards the Earth with a large velocity. Due to the 
Earth's magnetic field a net electrical current in the in- 
direction is induced with magnitude 



j{x, y, z,t)= / (v d q) e p e (x, y, z, t, h)dh , 



(3) 



where p~ e is the density of electrons and positrons, Eq. (2J. 
To take into account that the electrons and positrons 
(q = —1, +1 respectively) drift in opposite directions un- 
der the influence of the magnetic field, the average side- 
wards drift velocity is weighted with the charge, denoted 
as (v d q). 

The radius of curvature of orbits of the electrons with 
an energy e e = 7?nc 2 in the Earth's magnetic field is 
Rb = P"fmc/{eBE). A realistic magnitude of the mag- 
netic field {Be = 0.3 x 1CT 4 T) yields a curvature radius 
Rb = Pi x 50 m. The angular deflection is thus 9 = 
Lj Rb where L is the mean free path, i.e. the length over 
which the electrons scatter over a large angle due to mul- 
tiple soft scattering or a hard scattering. The transverse 
component of the velocity is vt = c6 — cLeB e / '(Pjmc) , 
assuming that sin 8 <C 1 or a transverse velocity much 
smaller than the longitudinal component. The drift ve- 
locity, being the average over the complete trajectory, is 
half this value 201. 



Vd 



c6/2 = 



cLcBe 
2Pjmc 



(4) 



The problem is thus now reduced to the calculation of 
the mean path length L. 

At high energies, e e > 10 MeV, the electron cross sec- 
tion is dominated by hard collisions and the mean path 
length is given by Lr ■ — Xq/ p a ir, where the electronic 
radiation length is X = 36.7 gem -2 and the density of 
air is p a i r — 10~ 3 gem -3 at sea level. The density is of 
course lower at higher altitudes. 
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FIG. 2: [color online] The drift velocity as function of shower 
age for different values of 70 in the multiple Compton scat- 
tering contribution, Eq. 0). The range in shower age corre- 
sponds to an elevation from, 10 km up to sea level. 



For smaller energies the above considerations do not 
apply since Coulomb scattering becomes the dominant 
scattering mechanism with a mean free path between 
successive collisions of 1 gcm~ 2 (which can be deter- 
mined directly from the corresponding cross sections). 
Since Coulomb scattering is strongly forward peaked, 
several collisions are necessary to randomize the momen- 
tum. As argued in Ref. [lj] the scattering angle due 
to multiple Coulomb scattering after traversing a dis- 
tance d, is S(9 2 } = d(^/^) 2 pair/ Xq with 70 = 40. When 
8(9 2 ) « 1, the original direction is lost, giving an effec- 
tive path length Lq = (7 /70) 2 X / p a i r - Over the whole 
energy range, the average path length L may be written 
as 1/L = 1/L C + 1/L R or 
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7- + 7o Pai 



(5) 



The drift velocity, obtained from Eqs. (|4l5[ ). now reads 

eB E X Q 



Vd 
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2 [3 7 2 + 7 mc fan- 



(0) 



keeping the energy dependence of the average path 
length. The average drift velocity is finally obtained by 
averaging over the energy distribution of electrons in a 
cosmic-ray air shower, using a parametrization given in 
Ref. [IH (see also Eq. (|A7[) ). Since this velocity is small, 
our assumption sin# <C 1 is indeed valid. 

The drift velocity depends rather strongly on the as- 
sumptions made in the estimate of multiple Coulomb 
scattering, as can be seen from Fig. [2] by choosing dif- 
ferent values for 70. The height dependence in (vaq) is 
due to the change in p a i r and due to the fact that the 
energy distribution of particles in the shower pancake de- 
pends on shower age [151 ]. In the present calculations we 
used a constant drift velocity, (vaq) = 0.04 c, equal to 
the value at the shower maximum for 70 =40. 



In Section IIIII the results obtained using the average 
current density is compared to the one obtained by ex- 
plicitly integrating the electric field generated by the par- 
ticles of different energies. The difference between the 
two appears to be mainly a normalization of the field 
strength, while the pulse shape is hardly affected. 



B. The vector Potential 

Given a current density, the vector potential can 
be obtained using the Lienard-Wiechert fields, 



1 



3» 



4tt£ J R(l-p-h) 



dh 



(7) 



ret 



for a source with an infinitcsimally small lateral exten- 
sion. We use the common notation where n is a unit 
vector pointing from the source to the observer and R is 
the distance, both evaluated at retarded time. Assum- 
ing that all particles move with the velocity of the shower 
front, the denominator in Eq. ([7]) can be rewritten to give 



V 



R(l - P a ■ h)\ ret = -(t- U) - n/3 s {h - C0 B t r ) 



y/(-cp a t + h) 2 + (1 - f3 2 n 2 )d 2 



(8) 



using Eq. (|B2|) for the retarded time. Eq. (jSJ) is written 
for a general medium with an index of refraction n how- 
ever all our calculations arc done in the limit n — 1. The 
distance between the observer and the point of impact of 
the core of the air shower is denoted by d, see Fig. Q] 

Since the current density has only an i-component, the 
vector potential will share this property, 



A x (t,d) =J dh 



,Pp{ h)ft{U) 
V 



(9) 



where T> is defined in Eq. ([5]) and where the current 
density is assumed to be parameterized according to 
Eq. ([3]) with J = (v^q) N e e/ A-ke c. We use SI units where 
= 1.44 x 10~ 9 [Vm] to get E in [V/m]. The ex- 
pression for the vector potential shown in Eq. ((§]) is the 
central equation in our derivation. 



1. Charge Conservation, Static Dipole 

At the point above the Earth's surface where the 
shower front passes, the charges arc being pulled apart by 
the Lorentz force. The air shower can thus be regarded 
as a 'zipper', pulling apart positive and negative charges 
at the point where it passes, leaving behind an electric 
dipole distributed along the path of the air shower. Since 
we have argued that the electric current, which is asso- 
ciated with the separating of the charges, is driving the 
electromagnetic pulse, we should also investigate the ef- 
fects of the created dipole. This dipole radiates because 
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it is not constant in time. To estimate its magnitude and 
the induced radiation field, we will assume that the pan- 
cake thickness is infinitely small. Please note that this 
dipole differs from the dipole mentioned in Q which is 
co-moving with the air shower. 

For definiteness, we temporarily assume that the 
charges are homogeneously distributed over a distance 
w in the x direction (this assumption will be relaxed 
at the end). For a shower front at an height z this 
corresponds to a line-charge density N e f t (—z/(3 s c)/w. 
Since the charges move sideways with a velocity (v^q), 
a charge Aq = ±At {vaq) N e f t {—z/ f3 s c)/w accumulates 
at x = ±w/2 after a time At. We will assume that this 
charge is at rest in the Earth's reference system and re- 
mains fixed at all later times while in reality it will slowly 
diffuse. Since the shower front progresses with a velocity 
(3 s c, vertical line-charge densities ±po(z) are created a 
distance w apart with 



Po(z) 



(v d g)N e f t (-z/f3 s c) 
c(3 s w 



(10) 



These charge densities at height z give a contribution to 
the zeroth component of the vector potential of magni- 
tude 

1 wp (z)x x -z 

" 47reo R 3 c[3 s R3 Jt{ p s c> ' [ ' 



d 2 



where we have introduced the distance R 
with d 2 = x 2 + y 2 and = ' z 2 + (x =F w/2) 2 + y 2 . 
From Eq. ([TTj) it is clear that the assumption of a homo- 
geneous line-charge density can be relaxed at this point. 
The scalar potential is obtained by integrating upward 
from the shower front over z, 



A°(t,d) = J^- 

CPs 



dz 



R 3 



(12) 



where at time t r the shower front has reached a height 
of zq = —(3 s ct r . The charges are now taken into account 
for the full development of the shower. It should, how- 
ever, be noted that gauge condition, d^A^, is not fulfilled 
since in the present simple model we have assumed that 
the charges forming the dipole moment are at rest in 
the Earth's system, while before they were moving with 
a vertical velocity @ s c. This sudden acceleration intro- 
duces an additional brcmsstrahlung contribution which 
is beyond the scope of the present work. 



2. Moving Dipole 

In the pancake, by virtue of the induced current, there 
will also be an induced electric dipole moving towards 
the Earth with the shower velocity. We will argue here 
that this dipole will not generate a contribution to the 
pulse in the limit used in this paper, n/3 s = 1. 



Due to the action of the Lorentz force the electrons and 
positrons will be displaced an average distance s. The 
contribution to the vector potential can now be written 
as 



A m A°(t,d) oc (—, — 

a?(l - 0ln 2 



8 1 



dxV 



V 3 



(13) 



using the same notation as introduced in Eq. (|11[) and 
calculate the derivative from Eq. ([8]). This contribution 
vanishes in the limit n(3 s = 1 and will be ignored in the 
following. 



C. The Electric Field 

The electric and magnetic fields can be derived from 
the vector potential in the usual way, 



E(t, d) = -d A(t, d) , 



(14) 



where we have ignored the zeroth component of the vec- 
tor potential (see the discussion at the end of this Sec- 
tion). 

Since the vector potential Eq. © has only a compo- 
nent in the x direction this will give rise to an electric 
field in the same direction. The emitted radiation is thus 
linearly polarized in the x direction, i.e. perpendicular to 
the shower axis and the magnetic field. 

The upper limit of the integral over h in Eq. (fT4f ex- 
tends up to infinity and we obtain 



E x {t,d) 



J / dh p p (h) 



dt jo 

d ft(tr) 



V 



dt V 



= -J 



J 



°d h m c0s d ^ 



V 



dh 



dh p p (h) 



d „ d \ ft{U) 
— + cp s — ] — 



dh) V 



c(3 s p p {h) 



Mr) 



V 



(15) 



h=0 



The second term can be rewritten as 

(d_ d_\ ft{t r ) 

\dt + Ps dh) V 

1 tdt T dt r \ df t (t r ) _ dft{t r ) 1 
V \ dt dh J dt r ~ dt r V 



(16) 



Using p p (h) = for h = 0, the expression for the electric 
field simplifies to 



E x (t,d) = -J [ 
Jo 



dh fjMpdp p (h) 



-J 



dh 



V ' dh 
p P {h) dftitr) 



V 



dt r 



(17) 
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where one should be careful in evaluating the integral 
because of the \j \Jcf5 s t — h pole in ^. To investigate the 
effect of this pole, associated with Cherenkov emission, 
we have explicitly studied the case (3 s n > 1 for which 
the pole in ^ lies inside the integration region. Only 
for unrealistically large values for the index of refraction 
Cherenkov radiation is emitted by the electric current 
density. For realistic values of n this effect is too small to 
distinguish. Since we see that our predicted pulse shapes 
for realistic values of n and n = 1 are identical we have 
limited ourselves to the latter. 



1. Limiting Case 

To obtain a simple estimate for the emitted radiation 
one may take the limit (3 S — 1 and n = 1 and ignore the 
thickness of the pancake, giving 



V = c/3 s t + (9(1 — (3 2 ) « ct 



and, for positive values of t, 
ct d 2 



ct r 



1 + (3 S 2c(3 s t 



+ 0(1 - ft) 



2ct 



(18) 



(19) 



which is large and negative since d 2 ^> c 2 t 2 . Interest- 
ing to note here is that the earlier part of the signal (t 
small and positive) contains the information of the earlier 
parts (at higher altitude) of the air-shower development 
(t r large and negative). The electric field can now be 
calculated, using J dh p v {h)h r* 0, J dh p p {h) = 1 



E x {t,d) = -J 



n 2 y /d 2 +c 2 2 t 2 r df t (t r ) 



Jft(t r ) 



cD 2 

cflt 
V 3 ' 



dtr 



(20) 



In the limit j3 s 
further to, 



E x (t,d) « j 



1 and n = 1 Eq. (f2"0"|) can be simplified 
c 2 t 2 A r. df t (u) 



rf 4 



dtr 



ft(U) 



(21) 



The limit t — > should be taken with care since this limit 
corresponds to large (and negative) retarded times, see 
Eq. (|19p . where the air shower may not even have started. 
As a result Eq. (|2Tj) produces a finite electric field at all 
times. For small distances to the shower core, d f» ct, the 
approximations made in deriving the expression for the 
retarded time, Eq. (|19p . are no longer valid and Eq. (|2"Tj) 
thus not applicable. 

From Eq. (|21|) it can be seen that the time structure 
of the pulse is, independent of the distance, given by a 
rather simple function of the longitudinal shower profile. 
If at distance d the peak of the pulse occurs at time to, at 
twice the distance, the signal peak occurs at a time 4 x t 
and the signal is four times as broad. From Eq. (|2"Tj) it 
can be seen that the peak value of the field, occurring 



at the same retarded time, has decreased by a factor 2 4 . 
It should be noted that the emitted radiation does not 
contain a rclativistic Lorcntz 7 factor and therefore does 
not depend on the exact velocity of the shower front, as 
long as it is close to c. The dependence on the energy 
distribution of the particles in the shower pancake is only 
indirectly through the dependence of the drift velocity. 

Eq. (|2~Tj) shows that there is a direct relation between 
shower profile and the pulse structure. The calculations 
presented in Section IIIII indicate that for a realistic case 
this feature is smeared out due to the effects of the finite 
extend of the pancake and the lateral distribution. 



2. Dipole Field 

To obtain an estimate of the effect of omitting the ze- 
roth component of the vector potential from our discus- 
sion the contribution due to the electric dipole, Eq. (| 1 2(1 . 
to the electric field is calculated as 



dA° 
dx 



= AE T = J 



nf t (-z/(3 s c)x' 2 



+J 



dz 



VR 3 
M-z//3 B c) z 2 +y 



z=z 
2 1 „,2 



2x 2 



(22) 



n ^ x . Following a similar approach we ob- 
tain for the other components 



using ^ 



J 



nf t (-z/(3 s c)xy 



+J 



VR 3 

dz 



z=z 
M-z/PsC) " 



3xy 



AE,. 



J 



0s 

nft(-z/p a c)xz 



R 5 



-J 



VR? 



dz 



z=z 

ft(-z/p s c) 3xz 



0s 



R 5 



(23) 



(24) 



It is important to note that in Eqs. (|22H24p the dis- 
tance R appears in the denominator in stead of T> as in 
Eq. (fl~7|) . The reason for this is that the electric dipole is 
at rest in the frame of the observer. Since R 3> T> for the 
cases of practical interest, the contribution of Eq. f22]) is 
much smaller than that of Eq. p7|) and thus can safely 
be ignored. 



D. Azimuthal Distribution of Radiation Pattern 

As remarked before, the emitted electric field due to 
the induced electric current, Eq. (|17|) . is linearly polar- 
ized in the x-direction. Its magnitude depends only on 
the distance to the shower core and, for a shower with a 
cylindrical symmetry, has a perfect azimuthal symmetry 
around the point of impact of the shower. This sym- 
metry is broken by the fact that, due to the drift veloc- 
ity which is induced by the magnetic field, the distribu- 
tion of charged particles in the shower is somewhat more 
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stretched in the ^-direction than in the y-dircction. In 
addition, the field of the induced dipole, Eq. (p?!?|) and 
Eq. (|23|) . docs not have an azimuthal symmetry. The 
symmetry breaking induced by these two effects is how- 
ever small and will not be considered further. 



III. RESULTS 
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Including a finite thickness of the pancake, using 
Eq. (fT?) with Eq. (|A5[) . results in the curves labeled 
'thickness' in Fig. [3] This has a very important effect 
on the pulse shape at 300 m, which can easily be un- 
derstood since a finie pancake thickness, L = cAt, intro- 
duces a 'smearing' effect for the pulse over a time At. 
Since for larger distances the pulse width is already siz- 
able (it increases roughly with the second power of the 
distance as follows from Eqs. (fT9]) and (f21) ) the smearing 
has only a minor effect at 700 m. 

Taking in to account the effects of the lateral spread of 
the particles in the shower, Eq. (|A6[) . in addition to the 
pancake thickness, results in drawn curve labeled 'lateral' 
in Fig. [3] The width of the pulse increases even further 
although the effect is relatively small. 

The effect of the static dipole field is shown in Fig. 21 
where it is compared with the pulse including the ef- 
fects of a finite pancake thickness, for two distances from 
the core. It should be noted that the dipole pulse has 
been multiplied by a factor 10 3 in order to be able to 
show it on the same scale. This clearly shows the asser- 
tion made earlier that the dipole response can be safely 
ignored. Apart from the small magnitude, also the as- 
sociated long wave length will make it undetectable in 
realistic experiments. 



FIG. 3: [color online] Effects of refinements of the shower 
structure on pulse shape at 300 m and 700 m from the shower 
core for a 10 17 eV shower. The line labeled 'appx' corresponds 
to the limiting case of Eq. H21\) while the other curves include 
the effects of the finite pancake thickness, the lateral extent 
of the shower and the energy distribution of electrons in the 
shower as discussed in the text. 

In Fig. [3] the calculated pulses as a function of time are 
shown at distances of 300 m and 700 m from the shower 
core for different levels of sophistication in the air-shower 
parametrization. The dotted curve, labeled 'appx', is the 
result of the most simple calculation using Eq. (f2Tj) where 
the longitudinal profile is given by Eq. (|A4[) . It has been 
verified that this result is indistinguishable from that ob- 
tained using the full expression, Eq. (|17p . in the limit of 
vanishing pancake thickness. In order to investigate the 
accurateness of this simple result as compared to a more 
realistic calculation we relax some of the approximations 
to see their effects. 

In arriving at Eq. (fT?) . the sideways drift velocity of 
the electrons and positrons has been averaged over their 
energy distribution. In doing so, the dependence of the 
denominator of Eq. ([TJ on electron energy has been ig- 
nored. To test the effect of this approximation we have 
instead used the full expression, see Eq. (|A8[) . resulting 
in the dash-dotted curves labeled 'velocity' in Fig.[3J still 
in the limit of vanishing pancake thickness. It shows that 
the effects of a spread in the energies of the electrons and 
positrons in the pancake affects mainly the magnitude of 
the pulse and hardly its time structure. At small dis- 
tances the decrease of the peak height due to this effect 
is stronger than at large distances. 
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FIG. 4: [color online] Comparing the pulse due to the dipole 
field with that of the current at two distances from the shower 



For a realistic shower one should expect a strong corre- 
lation between charged particle velocity and the distance 
behind the shower front (slower particles trailing further 
behind) . For this reason we have chosen not to mix the ef- 
fects of the velocity distribution and finite pancake thick- 
ness in the present work which is based on using simple 
parameterized showers. Results for a full Monte-Carlo 
simulation will be presented in a future work, which will 
also take into account the effects of an angular spread of 
the particles in the pancake. 

The effect of a finite pancake thickness can also be 
seen in the frequency response of the pulse as shown in 
Fig. O The frequency response is normalized such that 
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FIG. 5: [color online] Frequency response of the pulse at two 
distances from the shower core. In the calculations labeled 
'appx' and 'thickness' the thickness of the pancake is zero 
(Eq. \21\l ) and 10 ra respectively. 

j\E(t)\ 2 dt = 2 f\£(i>)\ 2 dv. At shorter distances d, the 
effect of finite thickness is to suppress the higher fre- 
quency components since the signal is only coherent for 
wave lengths larger than the typical size of the emitting 
system. The signal, in the limit where the thickness is 
ignored, does depend strongly on the distance from the 
shower core since the projected longitudinal extend en- 
ters, which equals to zero when viewing the shower head- 
on. Including a finite thickness reduces the dependence 
of the pulse shape on the distance from the core. 
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FIG. 6: [color online] Electric field strength (multiplied by d 2 ) 
at different distances from the shower core as function of time. 

In Fig.[S]the electric field is plotted as function of time 
for an observer at various distances d from the shower 
core. The primary energy is 10 eV and the calculation 



includes the effects of the pancake thickness only. The 
shower core hits the Earth's surface at t = 0. At large 
distances the pulse decreases in magnitude even faster 
than d 2 , as predicted by Eq. (|2"Tj) . At small distances im- 
portant deviations from the simple parametrization are 
observed. This is to a minor extent due to the fact that 
the approximations made to arrive at Eq. (|2ip are no 
longer valid, and mostly due to the effects of taking into 
account the finite thickness of the pancake which strongly 
influences the pulse shape at distances d < 500 m. In- 
cluding lateral extent of the shower will not greatly alter 
the picture. 
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FIG. 7: [color online] Fourier components of the electric field 
strength at different distances from the shower core for a 
10 17 e V shower. 
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FIG. 8: [color online] Same as Fig.^for a 10 19 eV shower. 

The frequency decomposition of the pulses shown in 
Fig. [6] are shown in Fig. [7] and those for an energy of 
10 19 eV in Fig. [S] At higher energies the shower maxi- 
mum is closer to the surface of the Earth making for a 
broader pulse. This is reflected in the frequency spec- 
trum by a peaking of the response at lower frequencies. 
At the same time the number of charged particles in the 
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shower is roughly proportional to the energy of the pri- 
mary particle. This in turn implies that the electric field 
is two orders of magnitude larger for a 10 19 eV induced 
shower than for 10 17 eV. 

The present results can be compared with those given 
in Refs. [l2|, [l3[. The basic features and magnitudes of 
the frequency responses are very similar. A difference is 
seen in the time structure of the peak. The pulse form 
obtained in this work has a distinct bi-polar structure, as 
has been discussed before, while that of Refs. [H,[ll| has 
a simple unipolar structure. The bi-polar structure of the 
pulse can also be understood from the fact that the vector 
potential is positive definite and vanishes for both small 
and large times. The electric field is the time derivative of 
this vector potential and crosses zero at the time when 
the vector potential reaches a maximum. For the first 
part of the pulse the two terms in Eq. (|21[) add construc- 
tively resulting it the large leading positive part. Since 
the vanishing vector potential at t = and t — oo does 
not depend on details of the shower profile (radial dis- 
tribution or velocity distribution) the predicted bi-polar 
shape, with a vanishing time-integral or zero-frequency 
component, can be regarded as a robust prediction. This 
prediction also follows from the work of ref. Q. The 
main difference between this work and ref. [|| lies in the 
fact that we have considered a more realistic shower pro- 
file and have presented a calculation directly in the time 
domain. The latter allowed us to show explicitly the re- 
lation between shower profile and pulse shape. 



IV. SUMMARY 

In this work a relatively simple macroscopic picture 
is presented for the emission of coherent electromagnetic 
radiation from an extensive air shower initiated by a high- 
energy cosmic ray. In this picture the radiation is emitted 
through the electromagnetic current which is induced by 
the Earth's magnetic field in the plasma at the front end 
of the shower. Under some simplifying assumptions a 
simple algebraic equation can be derived, Eq. (j!?T]) . which 
clearly shows which are the important aspects of the air 
shower that determine the electromagnetic pulse. 

It is shown that the time structure of the pulse directly 
reflects the longitudinal development of the number of 
electrons (and positrons) in the shower. The radio pulse, 
therefore, gives very similar information on the air shower 
as is obtained from air-fluorescence detection. Since the 
zero cross-over point of the pulse is related to the maxi- 
mum in the shower profile this implies that the peak in 
the pulse is related to the shower development well before 
the maximum. 

In this first paper we have restricted ourselves to a very 
simple geometry In a future publication emission from 
showers for a more general geometry will be investigated. 
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APPENDIX A: AIR-SHOWER 
PARAMETRIZATION 

The front of the shower moves with a velocity (3 s c in 
the — z direction, towards the Earth were the point 2 = 
is taken at the Earth's surface. At time t = the shower 
reaches the surface of the Earth. The shower thus exists 
at negative times only. 

The number of charged particles in the shower is pa- 
rameterized as function of time (t), distance from the 
shower front (h), and lateral distance (r = y 1 x 2 + y 2 ) 
where we consider here vertical showers only. For sim- 
plicity we assume that the different dependencies simply 
factorize 

p e (x,y,z,t,h) = N e ft(t)p KKG (x,y)p p (h)6(z + /3 s ct + h) . 

The distributions p p (h) and p NKG (r) [l6| are normalized 
such that their integrals equal to unity. The total number 
of charged particles at a specific time is thus given by 
N e ft(i). This simple paramctrization of the shower is 
sufficient to gain insight into the basic structure of the 
emitted electromagnetic pulse. 



1. Longitudinal Profile 

Following Ref. @, the longitudinal shower develop- 
ment can be parameterized using a shower age, 



s(X) 



3X/X 



X/Xo + 2X mx /X 



(Al) 



where .Xo=36.7 g/cm 2 is the electronic radiation length 
in air. The primary energy is denoted by E p and X the 
penetration depth in units of [gem ]. The parameter 
X mx is chosen such as to reproduce the positions of the 
shower maxima as have been determined from shower 
simulations fl~7j | . 



(A2) 



X mx = (840 + 70 log 10 {E p /10 20 eV) gem 



In the present calculations we model the atmosphere as 
X(z) = X(0)e- Cz (A3) 
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with X(Q) = 1000 gem 2 and where C is chosen such 
that X(ikm) = 630 genr 2 . Following Ref. @ the time 
dependence of the number of charged particles is param- 
eterized as 



NJ t (t) =N e e( x - x ™- 1 - 5X1 » s V Xo 



(A4) 



using Eq. (|A3j) with z = —/3 s t. The maximum number of 
charged particles is chosen as N e = 6 x (E p /10 10 eV) to 
agree with the number given in [17J for a 10 19 eV shower. 



Pancake Thickness 



In Ref. the measured arrival time distribution at 
a given radial distance is fitted with a r-probability dis- 
tribution function (T-pdf). Converted into a thickness of 
the shower front this can be written as Q 



p p {h) = h?e- 2h/L x (4/L 2 ) 



(A5) 



where the parameters /3=1 and L depend on shower age. 
At the shower maximum a reasonable choice for the pa- 
rameters is given by [3=1 and L = 10 m. In the present 
calculation we keep these fixed for the full development 
of the shower. The effects the curvature of the pancake 
has been ignored. 



3. Lateral Distribution 

The lateral particle density can be described with the 
NKG (Nishimura-Kamata-Greisen) [l6[ parametrization, 
which at the shower maximum (s = 1) reads 



. . 1 2.5 / r \- 1 / r 
Pnkg(2;,2/ = — — (H 



(A6) 



normalized such that 2ir / °° r drp NVLG (x,y) = 1. Here 
tm is the Molierc radius at the atmospheric height of the 
maximum derived from the atmospheric density as (l9| 

r Ai{h) = 9 p 6 ^ C "/ t ) ■ The atmospheric density at a height 
of 4 km corresponds to p a tm = 0.82 mgem , which in 
turn yields rjj ~ 117 m. 



taken to be 
from Ref. LL 



£;(e e ), where the parametrization is taken 



dpE(e e 



(A7) 



dloge e (e e + a)(e e + (i) s 

This parametrization is based on a detailed comparison 
with the results of Monte-Carlo calculations. The pa- 
rameters are taken as [l5[ a = 6.42522 — 1.53183s and 
= 168.168 - 42.1368s in units of McV where we have 
implemented this parametrization for shower age s = 1. 
The normalization constant Af is chosen to normalize the 
integral to unity. 

The energy distribution of the electrons and positrons 
can be included in the present calculation at two lev- 
els of sophistication. The simplest is to use it only in 
the calculation of the average drift velocity as discussed 
in Section III Al In a somewhat more sophisticated ap- 
proach the contribution to the vector potential, Eq. ([9]), 
is averaged resulting in 



.4' 



PE{(-e)vd{P) 

R + PP s t r c 



N e f t (t r )d\oge e 



(A8) 



where t r is negative and the shower is at height zq = 
—(3 s t r c. The denominator can be rewritten as R + 
(}(3 s t r c = U+zo(j3 s — f3). The shower front moves with ve- 
locity [3 s c while the charged particles move with a smaller 
velocity and thus must be trailing behind the front at a 
certain distance. The latter effect has not been included 
here. From the vector potential the electric field can be 
calculated in the usual way. 



APPENDIX B: RETARDED TIME 

At the retarded time t r the front of the shower is at a 
height z = —c[3 s t r . The travel time for a signal emitted 
at a distance h behind the shower front at z to reach the 
observer (at Earth at distance d from the point of impact 
on Earth) is 



t-t T = M = n^d 2 + (z + h) 2 /c , 



(Bl) 



where n is the index of refraction. The expression for the 
retarded time can thus be written as 



4. Energy Distribution 

In addition to the spatial distribution of the electrons 
in the shower, we have also taken their spread in energy 
into account. All the time we assume that the shower 
front is moving towards the Earth with velocity (3 S ~ 1. 
The distribution of electron energies e e in the shower is 



ct r 



ct - n 2 /3 s h - ny /(-c/3 s t + h) 2 + (1 - [3 2 n 2 )d 2 
1 - f3 2 n 2 



(B2) 

where in our convention the retarded time is negative for 
a shower above the Earth and t is assumed to be positive 
(the pulse, at a certain distance d, arrives only after the 
shower has hit the ground). 
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